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The  performance  of  polymer  electrolyte  membrane  fuel  cells  is  highly  influenced  by  the  water  content  in 
the  membrane.  To  prevent  the  membrane  from  drying,  several  researchers  have  proposed  extra  humidifi¬ 
cation  on  the  input  reactants.  But  in  some  applications,  the  extra  size  and  weight  of  the  humidifier  should 
be  avoided.  In  this  research  a  control  technique,  which  maintains  the  relative  humidity  on  saturated  con¬ 
ditions,  is  implemented  by  adjusting  the  air  stoichiometry;  the  effects  of  drying  of  membrane  and  flooding 
of  electrodes  are  considered,  as  well.  For  initial  analysis,  a  mathematical  model  reveals  the  relationship 
among  variables  that  can  be  difficult  to  monitor  in  a  real  machine.  Also  prediction  can  be  tested  optimizing 
time  and  resources.  For  instance,  the  effects  of  temperature  and  humidity  can  be  analyzed  separately.  For 
experimental  validation,  tests  in  a  fault  tolerant  fuel  cell  are  conducted. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Major  efforts  to  reduce  greenhouse  gas  emission  have  increased 
the  demand  for  pollution-free  energy  sources.  Fuel  cells  have 
attracted  great  attention  in  the  recent  years  as  a  promising  replace¬ 
ment  for  traditional  stationary  and  mobile  power  sources.  They 
are  characterized  by  high  power  density,  high  efficiency  and  low 
emissions. 

Significant  improvements  in  polymer  electrolyte  membrane 
fuel  cell  (PEMFC)  technology  have  been  achieved  over  the  past 
decade.  However,  the  advances  in  performance,  stability,  reliability, 
and  cost  for  today’s  fuel  cell  technology  are  not  enough  to  replace 
internal  combustion  engines.  A  number  of  fundamental  problems 
must  be  overcome  to  improve  their  performance  and  reduce  their 
cost. 

The  design,  control,  and  optimum  operation  of  fuel  cells  require 
an  understanding  of  their  dynamics  when  changes  in  current,  volt¬ 
age,  or  power  are  requested.  Also  a  control  system  is  needed  to 
ensure  that  humidity  and  temperature  are  within  prescribed  limits 
during  operation,  including  conditions  of  variable  loads. 

The  performance  of  a  PEMFC  is  highly  influenced  by  the  water 
content  in  the  membrane,  which  can  be  controlled  through  the  rel¬ 
ative  humidity.  In  general,  PEMFC  with  extra  humidification  of  the 
input  reactants  work  more  efficiently  [  1  ] ;  however,  in  some  applica¬ 
tions,  the  extra  size  and  weight  of  the  humidifier  should  be  avoided. 
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In  this  research  a  relative  humidity  control  technique  without  extra 
humidification  and  based  on  the  adjustment  of  the  stoichiometry  is 
implemented.  Initially,  the  control  technique  is  applied  on  a  PEMFC 
model  to  analyze  and  to  verify  predictions  about  the  system  evolu¬ 
tion.  After  that,  tests  are  conducted  for  experimental  validation  of 
the  simulation  results. 

This  paper  is  organized  as  follows.  In  Section  2,  the  basic  con¬ 
cepts  for  the  mathematical  model  of  a  PEMFC  are  introduced. 
Section  3  introduces  the  proposed  control  technique  and  presents 
simulation  tests.  Section  4  presents  tests  for  experimental  valida¬ 
tion  and  description  of  the  experimental  equipment.  In  Section  5 
conclusions  are  emphasized. 

2.  The  fuel  cell  model 

Several  mathematical  models  of  PEMFC  can  be  found  in  litera¬ 
ture  [1-4].  Basically  a  PEMFC  model  consists  of  an  electrochemical 
and  a  thermo-dynamical  part.  In  Ref.  [2],  an  electrochemical 
model  was  introduced;  in  this  case  the  polarization  curve  obtained 
with  this  model  is  compared  to  the  polarization  curve  of  the 
manufacturing  data  sheet  to  validate  the  model.  In  Ref.  [5],  the 
thermo-dynamical  part  of  the  model  is  introduced  and  the  electro¬ 
chemical  part  is  extended  to  consider  the  effects  of  different  types 
of  faults. 

2.1.  The  electrochemical  model 

The  electrochemical  model  permits  to  calculate  the  cell  voltage. 
The  output  voltage  VFC  of  a  single  cell  can  be  defined  as  the  result 
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Nomenclature 

A 

area  of  membrane  (cm2) 

^Nernst 

open  circuit  voltage 

FTFC 

fault  tolerant  fuel  cell 

he 

electrical  current  of  the  FC  (A) 

J  max 

maximum  density  of  current  (A  cm2) 

Jn 

fuel  crossover  (A cm2) 

£ 

thickness  of  membrane  (cm) 

nr 

number  of  cells 

PEMFC 

polymer  electrolyte  membrane  fuel  cell 

Ph2 

hydrogen  pressures  (atm.) 

Po2 

oxygen  pressures  (atm.) 

Psat 

saturated  vapor  pressure 

Pw 

partial  pressure  of  water  (atm). 

Qgen 

generated  heat  (W) 

RH 

relative  humidity 

Rc 

electrode  resistance  (£2) 

Pm 

membrane  resistance  (£2) 

7 

temperature 

Vact 

activation  voltage  drop 

Vcon 

concentration  voltage  drop 

VFc 

voltage  of  a  single  fuel  cell 

Vohmic 

ohmic  voltage  drop 

Vs 

stack  voltage 

Greek  letter 

X 

air  stoichiometric  relationship 

of  the  following  expression  [1  ]: 


Vpc  —  pNernst  —  Vact  —  ^ohmic  —  VCOn 


(i) 


^Nernst  represents  the  reversible  open  circuit  voltage: 


Ewernst  =  1.229-0.85  x  10“3(T  -  298.15)  +4.31 


x10-5T 


ln(PH2)+ 2  ln(Po2 


(2) 


where  Ph2  andP02  (atm)  are  the  hydrogen  and  oxygen  pressures, 
respectively,  and  T  (K)  is  the  operating  temperature. 

Vact  is  the  voltage  drop  due  to  the  activation  of  the  anode  and 
the  cathode: 


Vact  =  -[?!  +§2r  +  §3Tln(co2)  +  §4nn(/Fc)]  (3) 


Vcon  represents  the  voltage  drop  resulting  from  the  mass  trans¬ 
portation  effects,  which  affects  the  concentration  of  the  reacting 
gases: 

Vcon  =  -Bln(l-2_)  (6) 

where  B  (V)  is  a  constant  depending  on  the  type  of  cell,  Jmax  is 
the  maximum  electric  current  density,  and  J  is  the  electric  current 
density  produced  by  the  cell.  In  general,  J=Jout  +Jn  where  J0 ut  is 
the  real  output  electric  current  density,  and  Jn  represents  the  fuel 
crossover  and  internal  loss  current. 

Fig.  1  illustrates  the  polarization  curve  of  a  PEMFC  applying  the 
above  electrochemical  model. 

Considering  a  stack  composed  by  several  cells,  the  output  volt¬ 
age  can  be  assumed  to  be  Vs=nrVFC,  where  nr  is  the  number  of 
cells  composing  the  stack.  Flowever,  constructive  characteristic  of 
the  stack  such  as  flow  distribution  and  heat  transfer  should  be  taken 
[6-10]. 

2.2.  The  thermo-dynamical  model 

The  calculation  of  the  relative  humidity  and  the  operating  tem¬ 
perature  essentially  compose  the  thermo-dynamical  model. 


where  §,■  (i  =  1, . . .,  4)  are  specific  coefficients  for  each  type  of  cell,  JFc 
(A)  is  the  electric  current,  and  c02  (atm)  is  the  oxygen  concentration. 

Vohmic  is  the  ohmic  voltage  drop  associated  with  the  conduction 
of  protons  through  the  solid  electrolyte  and  electrons  through  the 
internal  electronic  resistance: 


V0hmic  =  +  Rc) 


(4) 


where  Rc  (£2)  is  the  resistance  to  electron  flow,  and  RM  (£2)  is  the 
resistance  to  proton  transfer  through  the  membrane: 


Rm  = 


Pm^ 


Pm  = 


181 .6  •  [l  +  0.03(/fc/A)  +  0.062(7/303  )2(/FC/A) 
[\/f  -  0.634  -  3(/fc/A)]  exp[4.18(7  -  303/7)] 


2.51 


(5) 


where  pm  (£2  cm)  is  the  specific  resistivity  of  membrane,  £  (cm)  is 
the  membrane  thickness,  A  (cm2)  is  the  active  area  of  the  mem¬ 
brane,  and  iff  is  a  specific  coefficient  for  each  type  of  membrane. 


2.2.7  Temperature 

The  variation  of  temperature  in  the  FC  is  obtained  with  the  fol¬ 
lowing  differential  equation  [2]: 


dT  AQ 
d t~MCs 


(7) 


where  M  (kg)  is  the  whole  stack  mass;  Cs  (JK_1  kg-1)  is  the  aver¬ 
age  specific  heat  coefficient  of  the  stack;  and  AQ  is  the  rate  of 
heat  variation  (i.e.  the  difference  between  heat  generated  and  heat 
removed).  Three  types  of  heat  removed  are  considered:  heat  by  the 
reaction  air  flowing  in  the  stack  (Qremi ),  by  the  refrigeration  system 
(Qrem2),  and  heat  exchanged  with  the  surroundings  (Qrem3)- 

The  rate  of  heat  generated  in  the  stack  is  calculated  from  the 
next  equation  [1]: 


Qgen  —  P0Ws 


(3) 


where  Pows  is  the  power  produced  by  the  stack. 
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2.2.2.  Relative  humidity 

The  proton  conductivity  is  directly  proportional  to  the  water 
content.  Thus,  there  must  be  sufficient  water  content  in  the  polymer 
electrolyte  membrane. 

In  this  research  extra  humidification  of  the  input  reactants  is 
avoided.  Then,  dry  hydrogen  comes  to  the  anode  and  air  at  envi¬ 
ronmental  conditions  comes  to  the  cathode. 

The  chemical  reactions  forms  water  at  the  cathode,  and  because 
the  membrane  electrolyte  is  very  thin,  water  would  diffuse  from 
the  cathode  side  to  the  anode  by  back-diffusion  during  the  opera¬ 
tion  of  the  cell.  On  the  other  hand,  the  hydrogen  protons  moving 
from  the  anode  to  the  cathode  pull  water  molecules  with  them  by 
electro-osmotic  drag.  Fortunately,  the  more  drag  the  more  water 
production  [11].  These  combined  effects  are  considered  as  the  net 
water  flux.  In  Refs.  [12-14],  the  net  water  flux  and  other  parameters 
are  considered  under  several  operational  conditions. 

The  water  production  would  keep  the  electrolyte  hydrated.  This 
level  of  hydration  can  be  measured  through  the  relative  humidity 
of  the  output  air.  In  practical  applications,  monitoring  the  rela¬ 
tive  humidity  of  the  output  air  is  a  very  simple  and  non-intrusive 
method  to  measure  the  relative  humidity  within  the  cell.  In  Ref.  [15] 
several  methods  for  monitoring  the  water  content  are  discussed, 
but  most  can  be  applied  only  in  laboratory  conditions. 

If  the  relative  humidity  is  too  low,  then  the  membrane  dries 
out  and  the  conductivity  decreases.  In  contrast,  a  relative  humidity 
too  high  produces  accumulation  of  liquid  water  on  the  electrodes, 
which  can  become  flooded  and  block  the  pores;  this  makes  gas 
diffusion  difficult  [14,16-19].  The  result  of  these  two  constraints  is 
a  fairly  narrow  range  of  operating  conditions. 

To  calculate  the  relative  humidity  of  the  output  air,  the  balance 
of  water  is  establishes: 


output  =  input  +  internal  generation 


or  in  terms  of  the  water  partial  pressure:  Pwout  =  ^win  +  Pwgen- 
And,  also  RH0ut  x  Psat.out  —  Pwout- 
then,  the  RH0Ut  is 


ou  pwin+ Pwge„ 

KHout  =  - p - 

1  sat_out 


(9) 


where  Pwin  is  the  water  partial  pressure  in  the  input  air;  Pwgen  *s  th e 
water  partial  pressure  generated  by  the  chemical  reaction;  Psa t_out 
is  the  saturated  vapor  pressure  in  the  output  air. 


Pwin  —  Psat_in^Hjn 

where  RHin  is  the  relative  humidity  of  the  input  air. 

PWgen  is  calculated  from  the  next  equation  [2]: 

_  42.1Pajr 

-  X- 0.188  ^10-) 


where  Pair  is  the  air  pressure  (atm)  and  A  is  the  air  stoichiometric 
relationship.  The  saturated  vapor  pressure  (kPa)  is  calculated  from 
the  following  equation: 


Psat  =  Ta 


exp 


((b/T)  +  c) 

10 


(11) 


If  T>  273.15  (I<),  then  a  =  -4.9283;  b  =  -6763.28;  c  =  54.22. 

Fig.  2  illustrates  the  variation  of  temperature  and  relative 
humidity  for  constant  air  stoichiometry  relationships,  A  =  2,  A  =  4, 
and  A  =  8. 

The  air  stoichiometry  (A)  is  calculated  according  to  the  following 
equation: 


A  = 


42.  IP-,; 


RHdesPSat_out  —  Pwin 


0.188 


(12) 


Fig.  2.  Temperature  and  relative  humidity  for  constant  stoichiometries  (A  =  2, 4,  and 
B). 

where  RFIdes  is  the  desired  relative  humidity  to  maintain  saturated 
condition,  normally  between  80  and  100%  [13]. 

To  prevent  the  membrane  from  drying,  several  researchers  have 
proposed  extra  humidification  on  the  input  reactants.  In  Ref.  [1], 
several  methods  of  extra  humidification  are  described.  In  general, 
PEMFC  with  extra  humidification  works  more  efficiently,  between 
20  and  40%  [12,20].  But  in  some  applications  (such  as  portable 
electronics),  the  extra  size  and  weight  of  the  humidifier  should 
be  avoided.  In  application  of  middle  power  PEMFC  (about  5  kW) 
the  humidification  sections  accounts  for  about  20%  of  the  total 
volume  and  weight  [12].  And  most  likely,  the  smaller  plant  the  big¬ 
ger  humidifier.  Also,  some  researches  have  reported  other  forms  to 
avoid  extra  humidification  [13,21,22]. 

Nevertheless,  the  operating  conditions  of  the  PEMFC  are  more 
sensitive  to  variation  of  the  input  air  temperature  than  input  rel¬ 
ative  humidity.  The  effects  of  temperature  and  relative  humidity 
of  the  input  air  can  be  analyzed  separately  thanks  to  the  PEMFC 
model.  Fig.  3  illustrates  the  variation  on  the  output  RFI  produced 
by  the  adjustment  on  the  input  RH  when  input  temperature  is  con¬ 
stant.  In  Ref.  [13]  a  similar  analysis  is  made  considering  RHin,  RH0Ut, 
and  A. 

Another  method  to  confirm  the  reduced  effect  of  the  input 
RH  on  the  operating  conditions  of  the  PEMFC  is  to  verify  the 
limit  operating  temperature.  The  limit  operating  temperature  is 
the  temperature  in  which  the  PEMFC  can  operate  preserving  a 
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minimum  recommended  stoichiometry  (in  this  case,  A  =  2)  and  a 
recommended  output  RH  (in  this  case,  RH0Ut  on  saturated  con¬ 
ditions).  Fig.  4  illustrates  the  variation  of  the  limit  operating 
temperature  obtained  from  combining  Eqs.  (11)  and  (12)  when 
RHin  is  modified  and  the  temperature  of  the  input  air  is  constant 
(Tin-  25  °C). 

In  contrast,  Fig.  5  illustrates  the  variation  of  the  limit  operating 
temperature  when  Tin  is  modified  and  RHin  is  constant  (RFIin  =  50%). 
Actually,  the  water  partial  pressure  of  the  air  at  T=60°C  and 
RHinput  =  50%  is  more  than  three  times  higher  compared  to  this  at 
T=  25  °C  and  RFIinput  =  100%. 

To  prevent  flooding  problems,  a  common  method  for  remov¬ 
ing  excess  water  inside  the  PEMFC  is  using  the  reaction  air  flowing 
through  it.  Air  would  be  blown  over  the  cathode,  and  apart  from 
supplying  the  necessary  oxygen  it  would  dry  out  any  excess  water. 
The  correct  variation  of  the  stoichiometry  A  would  maintain  the  RFI 
proximal  to  saturated  conditions. 

The  air  stoichiometry  influences  both  the  availability  of  oxygen 
as  well  as  the  humidity  of  the  membrane.  A  low  stoichiometry  limits 
the  availability  of  oxygen  because  the  air  is  depleted  of  oxygen  when 
it  reaches  the  end  of  the  airflow  channels.  In  Ref.  [23],  the  effects  on 


Fig.  5.  Limit  operating  temperature  vs.  rinpUt. 


the  PEMFC  performance  with  different  level  of  fuel  utilization  and 
air  stoichometry  utilization  were  tested.  In  general  the  maximum 
efficiency  occurs  at  about  80%  of  fuel  utilization  (H2)  and  25%  of  air 
utilization. 

The  flooding  of  electrodes  makes  gas  diffusion  difficult  and 
affects  the  performance  of  the  PEMFC.  This  effect  changes  the  resis¬ 
tance  of  electrodes  (Rc),  and  is  simulated  applying  the  following 
equation  [5]: 


^c(k)  —  ^c(0) 


/  Wacum(/<) 

y  consti 


0.8 


(13) 


where  RC(0)  is  the  electrode  resistance  at  the  initial  state  (normal 
condition),  consti  is  a  constant  defining  when  the  electrodes  are 
led  to  flooding,  and  Wacum(fc)  is  the  volume  of  water  accumulated  at 
instant  k,  WacUm  =  mH2o(RH  -  100%).  The  rate  of  water  production 
(kgs-1)  is  calculated  from  the  following  equation  [1]: 

mH2o  =  9.34  x  10_8/FCnr  (14) 


This  effect  also  can  reduce  the  density  of  current;  this  reduction  is 
simulated  by  the  following  equation  [5]: 


7. 


max(k)  : 


J  max(0) 


(Wacum(fc) /consti 


A.2 


(15) 


where  7max(o)  is  the  value  of  the  maximum  electric  current  density 
at  the  initial  state  (normal  condition). 

On  the  other  hand,  when  temperature  is  higher  than  the  limit 
operating  temperature,  the  reacting  air  has  a  drying  effect  and 
reduces  the  RH.  A  low  RH  can  produce  a  catastrophic  effect  on  the 
polymer  electrolyte  membrane,  which  not  only  totally  relies  upon 
high  water  content,  but  also  is  very  thin  (and  thus  prone  to  rapid 
drying  out).  The  drying  of  the  membrane  changes  the  resistance  to 
proton  flow  (RM)-  is  affected  by  the  adjustment  of  (see  Eq. 
(5));  its  variation  is  simulated  according  to  the  following  equation 
[5]: 


^(fc)  = 


_ fm _ 

(const2/RHout(k))112 


(16) 


where  t/^o)  is  the  value  at  saturated  condition,  RH0Ut(k)  is  the  relative 
humidity  of  the  output  air  at  instant  k,  and  const2  is  a  constant 
defining  when  the  membrane  is  led  to  drying. 

Table  1  presents  the  parameters  of  the  PEMFC  in  normal  condi¬ 
tions. 

In  Ref.  [24],  a  multi-parametric  sensitivity  analysis  is  performed 
to  define  the  importance  of  the  accuracy  of  each  parameter.  The 
accuracy  was  analyzed  in  normal  conditions,  considering  variations 
around  ±10%  of  their  normal  values. 


Table  1 

Parameters  of  the  FC 


Parameter 

Value 

nr 

4 

A 

62.5  cm2 

i 

25  fxm 

po2 

0.2095  atm 

Pu2 

1.47628  atm 

Rm 

0.003  a 

B 

0.015  V 

£i 

-0.948 

£2 

0.00286  +  0.0002  \nA  +  (4.3  x  10-5)lncH2 

£3 

7.22  x  10-5 

£4 

-1.06153  x  10-4 

f(0) 

23.0 

Jn 

0.022  A  cm2 

J max 

0.672  A  cm2 
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Fig.  6.  Effect  of  the  relative  humidity  on  the  FC  internal  resistances. 


Also,  in  Refs.  [3,17]  the  variation  of  the  resistances  was  associ¬ 
ated  with  fault  detection  of  flooding  and  drying. 

Fig.  6  illustrates  the  influence  of  the  relative  humidity  on  the 
resistances  according  to  Eqs.  (15)  and  (16). 


Fig.  8.  The  FTFC. 


3.  The  fuel  cell  control  technique 

The  proposed  FC  system  adjusts  the  air  stoichiometry  to  main¬ 
tain  humidity  constantly  on  a  saturated  condition.  Flowever,  in  a 
real  machine,  there  is  not  direct  control  on  the  stoichiometry,  but 
there  is  direct  control  on  the  airflow. 


time  (min)  time  (min)  time  (min) 


time  (min)  time  (min)  time  (m'n) 


time  (min) 


Fig.  7.  Operation  of  the  PEMFC. 


10  20  30  40  50  60 

time  (min) 
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Fig.  9.  Monitoring  of  the  FTFC  applying  LabView. 


Applying  Eq.  (17)  we  calculate  the  reaction  airflow  (Ls-1 )  [1  ]. 

flow  =  3.0238  x  10_4/FC_adjnrA  (17) 

The  electric  current  (/FC)  has  been  adjusted  for  considering  internal 
loss  current  and  fuel  crossover  (Jn). 


In  general,  JFC_adj  =/Fc-reai  +  VL 

/pc-reai  depends  on  the  requested  load,  and  nr  is  the  number  of 
cells  (i.e.  a  constructive  parameter).  Therefore,  the  adjustment  of 
the  airflow  is  performed  through  the  adjustment  of  A  according  to 
Eq.  (12). 


Fig.  10.  Controller  block  diagram  in  LabView. 
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In  the  tests,  air  compressor  and  extra-humidification  system  are 
avoided.  Then  the  air  is  considered  with  pressure  of  1  atm  and  at 
environmental  temperature  (25  °C)  and  humidity  (50%).  Also  the 
PEMFC  system  operates  below  the  limit  operating  temperature.  And 
the  hydrogen  supply  is  dead-end. 

Before  experimental  validation,  analysis  on  the  PEMFC  model 
was  performed.  The  model  is  based  on  the  equations  of  previous 
section  and  the  parameters  presented  in  Table  1.  In  Ref.  [5]  this 
model  was  applied  to  verify  the  effects  of  different  types  of  faults, 
including  drying  and  flooding. 

Fig.  7  illustrates  the  evolution  of  the  main  PEMFC  variables  as  a 
function  of  time.  The  variables  are:  voltagestack  (V),  electric  current 
JFC  (A),  temperature  (°C),  volume  of  airflow  (Ls-1 ),  generated  heat 
(W),  stoichiometry  A,  power  (W),  and  relative  humidity.  Initially  in 
this  simulation,  the  PEMFC  supports  a  constant-load  demand;  thus, 
the  voltage  and  current  should  adjust  by  themselves  to  supply  this 
demand  (i.e.  the  output  power  would  be  constant).  And  the  control 
system  adjusts  the  airflow  to  maintain  the  humidity  on  the  desired 
value. 

The  simulation  begins  with  the  PEMFC  system  in  stand-by  (i.e. 
without  load,  and  at  environmental  temperature,  approximately 
25  °C).  After  the  load  requirement  is  done,  the  electrical  equi¬ 
librium  is  reached  in  less  than  3  s  (such  as,  the  equilibrium  of 
voltage  and  current).  On  the  other  hand,  the  temperature  increases 
slowly  as  a  consequence  of  a  high  inertia  of  the  thermo-dynamical 
state.  At  t  =  30min,  the  thermo-dynamical  state  is  almost  stable, 
then  step-variations  of  load  are  preformed  at  t  =  30  and  t  =  45min 
to  analyze  the  transient  response  (variation  of  20  and  50%  of 
load,  respectively);  in  these  cases  voltage  and  electric  current 
are  self-adjusted  for  the  requested  load,  the  airflow  is  adjusted 
by  the  control  system,  note  that  the  RH0Ut  remains  constant 
(85%). 

4.  Experimental  validation 

Key  issues  modeling  PEM  fuel  cell  systems  are  still  developing, 
those  include:  lack  of  measurement  techniques,  especially  real  time 
(in  situ)  and  non-intrusive  techniques.  More  work  is  required  in 


the  areas  of  modeling,  measurement  methods  and  fuel  cell  design 
optimization  [23]. 

A  fault  tolerant  fuel  cell  (FTFC)  was  constructed  at  the  PSERC 
laboratory  [25,26].  The  control  system,  the  sensor  system,  and  the 
power  system  compose  the  FTFC.  The  control  system  performs  the 
adjustment  of  the  airflow  blower  and  the  refrigeration  blower.  The 
sensor  system  implements  the  monitoring  of  voltage  (Vs)»  elec¬ 
tric  current  (JFC),  temperatures  (Tout  and  Tin),  and  relative  humidity 
(RHout  and  RHin).  The  power  system  is  composed  by  one  AvistaL- 
abs  cartridge  containing  four  polymer  electrolyte  membranes.  The 
software  LabView  executes  the  control  technique.  The  same  Lab- 
View  is  applied  for  monitoring  the  variables  and  for  controlling  the 
speed  of  the  blowers.  The  airflow  is  supplied  by  a  brushless  DC  axial 
blower,  nominal  feed  12  VDC,  maximum  airflow  3  CFM  (cubic  feet 
per  minute),  operating  voltage  5-13.8  VDC,  maximum  power  0.7  W, 
and  maximum  speed  4000  rpm.  The  hydrogen  supply  is  dead-end 
with  periodic  manual  purging  to  release  impurities. 

Fig.  8  illustrates  the  FTFC,  the  load,  and  a  desktop  computer 
with  the  software  LabView  executing  the  monitoring  process.  The 
constructive  parameters  of  the  FTFC  are  presented  in  Table  1 . 

Fig.  9  illustrates  the  graphic  user  interface  of  the  panel  control 
applying  the  LabView.  This  panel  control  presents  several  variables 
of  the  FTFC  such  as  output  voltage  (Vs),  electric  current  (JFC),  tem¬ 
perature,  relative  humidity  (RH0Ut),  and  airflow  volume.  In  Fig.  9 
the  FTFC  begins  to  operate  from  stand-by  and  powering  a  constant 
load  demand. 

The  control  technique  was  implemented  applying  LabView; 
then  the  equations  considered  on  previous  sections  should  be 
translated  to  the  LabView  syntax.  Fig.  10  illustrates  the  block 
diagram  of  the  controller  applying  LabView.  The  inputs  of  the  con¬ 
troller  are:  Jn,  area,  IF c_real,  Tin,  RHin,  voltage,  and  Tout.  Inside  the 
^Air  block,  Eqs.  (12)  and  (17)  are  applied  to  calculate  the  air  reac¬ 
tion  flow.  Inside  the  PFAN  block  (power  of  the  refrigeration  fan)  an 
empirical  equation  is  applied  to  calculate  the  power  of  the  refrig¬ 
eration  subsystem  (Pfan-Ref)-  ^fan-Ref  depends  on  the  voltage,  /Fc-adj 
and  Tout* 

Fig.  11  illustrates  the  inside  of  the  FAir  block,  basically  in  this 
block  Eq.  (17)  is  calculated  and  the  power  of  the  airflow  blower  is 


Fig.  11.  FAir  block. 
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Fig.  12.  Evolution  of  the  FTFC  variables  under  step  load  requirement. 


limited  by  the  maximum  air  volume,  3  CFM  (cubic  feet  per  minute). 
Eq.  (12)  is  calculated  inside  the  A,control  block,  in  a  similar  method. 

Before  recording  experimental  test,  the  thermodynamic  steady 
state  was  established.  In  Fig.  9  the  evolution  of  the  FTFC  variables 
to  reach  the  thermodynamic  equilibrium  was  illustrated.  After  the 
equilibrium,  the  experimental  test  considering  a  step  variation  of 
load  was  conducted.  In  this  test  the  load  was  reduced  by  50%  at 
t  =  30  min,  and  reestablished  at  t= 45  min.  Fig.  12  illustrates  the  evo¬ 


lution  of  output  voltage  (Vs)  and  electric  current  (JFC),  and  also  how 
the  controller  maintains  a  constant  relative  humidity  (RH0Ut). 

5.  Conclusions 

The  relative  humidity  is  a  crucial  variable  to  improve  the  perfor¬ 
mance  and  to  avoid  damage  to  the  membranes  of  PEMFC.  A  control 
technique,  which  maintains  a  satisfactory  relative  humidity,  was 
implemented.  This  technique  regulates  the  airflow  adjusting  the 
air  stoichiometry  and  therefore  maintains  the  relative  humidity  in 
a  recommended  value. 

A  PEMFC  model  was  applied  to  analyze  the  evolution  and  to 
verify  the  dependence  among  the  variables.  From  the  mathemat¬ 
ical  model,  the  evolution  of  some  variables  that  can  be  difficult 
to  monitor  in  a  real  machine  are  observed  (such  as,  generated 
heat,  stoichiometry,  portion  of  non-evaporated  water,  humidity 
over  100%,  etc.).  Also  tests  that  can  imply  in  permanent  damage  to 
the  equipment  can  be  avoided  ( such  as,  tests  in  very  dry  conditions). 

In  addition,  predictions  about  the  evolution  of  those  variables 
can  be  tested,  optimizing  time  and  resources.  The  analysis  in  the 
PEMFC  model  shows  that  the  relative  humidity  control  strategy  is 
stable  and  consistent  under  different  operational  conditions. 

Experimental  tests  were  conducted  to  verify  the  reliability  of 
the  proposed  technique.  The  technique  proved  to  be  easy  for  imple¬ 
mentation  in  the  FTFC  equipment. 
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